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ABSTRACT 

Aims. Our goal is to investigate the properties of the circumstellar disks around intermediate mass stars to determine their occurrence, 
lifetime and evolution. 

Methods. We completed a search for circumstellar disks around Herbig Be stars using the NRAO Very Large Array (VLA) and the 
IRAM Plateau de Bure (PdB) interferometers. Thus far, we have observed 6 objects with 4 successful detections. The results towards 
3 of these stars (R Mon, MWC 1080, MWC 137) were presented elsewhere. We present our new VLA and PdBI data for the three 
objects MWC 297, Z CMa, and LKHa 215. We constructed the SED from near-IR to centimeter wavelengths by adding our millimeter 
and centimeter data to the available data at other wavelengths, mainly Spitzer images. The entire SED was fitted using a disk+envelope 
model. In addition, we compiled all the disk millimeter observations in the literature and completed a statistical analysis of all the 
data. 

Results. We show that the disk mass is usually only a small percentage (less than 10%) of the mass of the entire envelope in HBe stars. 
For the disks, there are large source-to-source variations. Two disks in our sample, R Mon and Z CMa, have similar sizes and masses 
to those found in T Tauri and Herbig Ae stars. The disks around MWC 1080 and MWC 297 are, however, smaller (r oi „<100 AU). 
We did not detect the disks towards MWC 137 and LkHo- 215 at millimeter wavelengths, which limits the mass and the size of the 
possible circumstellar disks. 

Conclusions. A comparison between our data and previous results for T Tauri and Herbig Ae stars indicates that although massive 
disks (~0. 1 M Q ) are found in young objects (~10 4 yr), the masses of the disks around Herbig Be stars are usually 5-10 times lower 
than those around lower mass stars. We propose that disk photoevaporation is responsible for this behavior. In Herbig Be stars, the 
UV radiation disperses the gas in the outer disk on a timescale of a few 10 3 yr. Once the outer part of the disk has vanished, the entire 
gaseous disk is photoevaporated on a very short timescale (~10 5 yr) and only a small, dusty disk consisting of large grains remains. 

Key words, stars: individual: MWC 137, LkHa 215, R Mon, Z CMa, MWC 297, MWC 1080-stars: formation-stars :pre-main- 
sequence - circumstellar matter - planetary systems: formation - planetary systems: protoplanetary disks 
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■ 1. Introduction excesses and Group II sources with modest infrared excesses. 

' Group I sources have Chiang & Goldreich-like flar ed disks and 
Herbig Ae/Be (HAeBe) stars are intermediate-mass (M~2- Group n sources host flatter> se if. sh adowed disks. LeinertetalJ 
8 Mo) pre-main-sequence objects. Since these objects share (2004) confirmed this hypothes i s on the basis of 10 /on inter- 
many characteristics with high-mass stars (clustering, PDRs), feromet data while most HAe stars bel to G I an d 
but are far closer to us and less embedded, the detection of err- m surrounded b disks similar to those in T Tauri (TT) stars , 
cumstellar disks around these stars is crucial to the understand- Hgrbi Be (Hfie) stars arg mosd found in G n and seem 
ing of the massive-star-formation process. Herbig Ae (HAe) tQ be surrounded b disks with a flatter geometry dAcke et al.1 
stars are the precursors of Vega-type systems: determining the |2QQ5[) 
frequency and timescales of the disks around HAe stars is there- 
fore also important for planet formation studies. Although the existence of different types of disks is well ac- 
Substantial theoretical and observational efforts have been cepted, whether this difference is due to the different character- 
carried out recently for the under standing of the disk occurrence istics of the central star (hotter stars produce flatter disks), or 
and evolution in HAeBe stars dMeeus et alj l2001t IVink et al.1 to different characteristics of the circumstellar matter (different 
l2002t iMillan-Gabet et al.ll200H lAcke et al.ll2005l). On the basis grain opacities), or to an evolutionary link between flared, self- 
of the different mid-infrared excesses. iMeeus et al.1 d200ll) pro- shadowed, and debris disks is not established. If the evolution- 
posed that disks around Herbig Ae/Be stars can be classified into ary scenario stands, the different geometry between disks in HAe 
two Groups that they interpreted in terms of different geome- and HBe stars is the consequence of a more rapid grain growth 
tries: Group I sources with strong mid-infrared (20-100 //m) flux and a shorter dissipation timescale in the latter. The grain growth 
causes the optical depth of the disk to decrease and allows the 
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photo evaporate the disk external layers (iDullemond & Dom inik 
|2004|) . Disks around HBe would flatten and lose a large fraction 
of their mass before the pre-main-sequence phase (<0.1 Myr). 

One problem in discerning between the different scenarios is 
that most studies are based on the optical-NIR and mid-IR obser- 
vations. Optical-NIR and mid-IR observations probe the region 
within a few AU around the star. In addition, since the disk is 
optically thick at these wavelengths, they only provide informa- 
tion about the disk surface and cannot be used to derive the disk 
mass. An obvious step towards the formation of planetesimals 
in disks is coagulation and growth of the su bmicron sized grain s 
accreted from the proto-stellar envelope (Domi nik et al.l l2007). 
Although the end results are plainly visible in our own planetary 
system, the details and mechanisms of dust evolution are not un- 
derstood. The main problem is that the cooler, large grains do 
not emit significantly at infrared wavelengths. Only by observ- 
ing the submm-mm range of the spectrum are we able to deter- 
mine the sizes and properties of these large grains. Millimeter 
observations must be used to trace this more evolved grain com- 
ponent. Thus, imaging disks at mm wavelengths is required to 
trace the outer part of the disk and determine its mass and grain 
properties. 

We have carried out a search for circumstellar disks around 
HBe stars using the NRAO Very Large Array (VLA) and the 
IRAM Plateau de Bure (PdB) interferometers to investigate the 
properties of the circumstellar disks around intermediate-mass 
stars and to determine eventually their occurence, lifetime, and 
evolution. Some re suits have already been published el sewhere 
dFuente et alj2003l l2006t [Alonso-Albi et al.l2007ll2008l) . In this 
Paper, we present new observational data (LkHa 215, Z CMa 
and MWC 297) that complete our survey and we provide a com- 
prehensive analysis of all the data obtained thus far. 



2. Observations 

2.1. Selected sample 

The list of observed sources is shown in Table Q] Our sample 
was built according to the following criteria: Firstly, we included 
most well-studied Herbig stars with spectral type close to BO 
in the northern hemisphere (R M on, MWC 1080, M WC 137, 
MWC 297). In a pioneering work. iNatta et all (12000) compiled 
the existing 1.3mm/2.7mm interferometric disk detections in 
HAeBe stars. By this time, only one disk had been detected 
in stars with spectral type earlier than B7-8. They proposed 
that this lack of detections could be due to a shorter dissipa- 
tion timescale for the circumstellar disk in stars of spectral type 
earlier than B7. If the disk mass decreases with stellar spectral 
type, one expects to measure the minimum disk masses towards 
the most massive B0 stars. Secondly, we included some objects 
with spectral types between B7 and B8, where the transition be- 
tween the TT-like and massive star disks, if it exists, is expected 
to occur. Additional criteria were the closeness to the Sun and 
previous detections with single-dish telescopes. 

2.2. Millimeter and centimeter observations 

Thus far, we have observed 6 HBe stars with 4 success- 
ful detections. The VLA and PdBI observations towards R 
Mon, MWC 1080 and M WC 137 were already reported by 
iFuente et al.l ([2003, 2006). In this Paper, we present new VLA 
observations towards MWC 297 and Z CMa, and new PdBI 
observations towards LkHa 215, Z CMa and MWC 297. A 



Table 1. Sample of observed Herbig Be stars 



Name 


R.A. 


Dec. 


Dist. 


Mass 


Sp. 




(J2000) 


(J2000) 


(pc) 


(M Q ) 


type 


MWC 137 


06:18:45.50 


15:16:52.4 


1300 


14 


B0 


R Mon 


06:39:09.95 


08:44:10.7 


800 


8 


B0 


MWC 1080 


23:17:25.57 


60:50:43.3 


1000 


10 


B0 


MWC 297 


18:27:39.53 


-03:49:52.0 


250 


9 


B1.5 


LKHa 215 


06:32:41.79 


10:09:33.6 


800 


7 


B7.5 


ZCMa 


07:03:43.16 


-11:33:06.2 


930 


12 


B8 



summary of these observations, wavelength, beam, and date, is 
shown in Table [2] 

The data were calibrated, mapped, and analyzed in the 
GILDAS software package. Continuum images at mm wave- 
lengths have been produced by averaging the channels free of 
line emission after careful visual inspection. Natural weightings 
were applied to the measured visibilities producing the beams 
shown in Table [2] Since the disks are expected to be unresolved 
by our interferometric observations, we used the peak flux to 
construct the SED. The measured peak fluxes are shown in Table 
3. The rms of the images are shown in Tables A.l, A.3, A.5, A.7, 
A.9andA.ll. 

2.3. Mid-IR observations 

Four of the 6 sources in our sample have been observed at 24 /mi 
either with MIPS or IRS onboard the Spitzer Space Telescope, 
namely MWC 1080, MWC 297, R Mon, and LkHa 215. We 
retrieved the post-basic-calibrated data (pbcd) for imagery and 
basic-calibrated data for spectroscopy (bed). The pbcd files were 
used without additional processing. The bed spectra l mapping 
files were processed with the CUBISM software dSmith et al.l 
120071) . 

3. Results 

3.1. Interferometric data 

The results of our observations are shown in Table [3] We de- 
tected emission at 1.3mm and/or 2.7mm towards 5 out of the 6 
stars studied. Emission at centimeter (cm) wavelengths using the 
VLA array was detected towards Z CMa, MWC 297, MWC 137, 
and MWC 1080. In these massive stars, a significant fraction of 
the 1.3mm flux can be attributed to free-free emission instead 
of dust, thermal emission arising in the circumstellar disk. To 
estimate this contribution, we fitted the free-free spectral index 
at cm wavelengths and extrapolated the emission to millimeter 
(mm) wavelengths. In the case of MWC 137, all the flux seems 
to be due to free-free emission without any significant excess at 
1.3mm and 2.7mm. In the other cases, there is some flux excess 
at mm wavelengths that can be interpreted as arising in a cir- 
cumstellar disk. We show in Table [3] the spectral index between 
1.3mm and 2.7mm interferometric observations. Two values are 
presented: the first represents the total observed fluxes and the 
second with the free-free emission substracted, i.e., the spectral 
index of the thermal, dust emission. The 1.3mm/2.7mm spec- 
tral index is quite low (<2.2) for all sources. As discussed in 
the following Sections, this low spectral index constitutes obser- 
vational evidence of either optically thick dust emission or the 
presence of large grains in the disk midplane (see e.g. Alonso- 
Albi et al. 2008). As we discuss in Sect. 5.7, the second possibil- 
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Table 5. Grain properties in the disks 1 . 
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Name 


A (mm) 


Beam size 


Date 


Telescope 


MWr 7Q7 




1.1" x 0.4" 


Fph 7006 


PdBI 




~> 6 


1 4" y 0" 


Fph 7006 




T KH<t 7 1 5 


1 .3 


3.3" x 1.8" 


Anril 7007 


PdBI 




2.6 


6.8" x 3.7" 


April 2002 




ZCMa 


1.3 


2.3" x 0.9" 


Jan/March 2006 


PdBI 




2.6 


4.1" x 2.0" 


Jan/March 2006 




MWC 297 


7 


1.9" x 1.6" 


Dec 2005 


VLA 




13 


3.9" x 3.1" 


Dec 2005 




ZCMa 


7 


1.4" x 1.4" 


Oct/Dec 2005 


VLA 




13 


3.5" x 2.7" 


Oct/Dec 2005 






36 


9.2" x 6.6" 


Oct/Dec 2005 





Name 



Midplane 
grains 



Surface 
grains 



ity is the most plausible in the circumstellar disks around R Mon 
and Z CMa. 

In Table [3] we also show the ratio of our interferometric 
fluxes to previous single-dish data. In most cases, the flux de- 
tected by interferometric techniques is far lower (<10%) than 
that detected by single-dish observations. This implies that the 
envelope emission dominates the single-dish flux at FIR and mm 
wavelengths in most HBe stars. For this reason, interferometric 
measurements are essential to estimating the true disk emission. 
Since interferometric observations are unavailable at all wave- 
lengths, a disk+envelope model is required to fit the observed 
SED properly. 



3.2. Observed SED 

We reconstructed the entire SED towards our sample by com- 
pleting our VLA, PdBI, and Spitzer photometry with the data 
available in the literature. At optical and NIR wavelengths, ob- 
servations are heavily affected by extinction. To correct for this 
reddening, we substracted the standard B-V color correspon ding 
to the spectral type of each star, as given by iJohnsonl d 19661) . to 
the observed B-V color. This color excess E was used to calcu- 
late the visual extinction A v = R y E assuming a ratio of total-to- 
selective absorption of R v = 3. 1 dCardelh etaljl989h . With these 
two parameters, we estimate the extinction corresponding to a 
particular obs e rvatio n using the parameterized extinction law of 
ICardelli et alJ (fl989h . which can be applied to wavelengths be- 
tween 0.12 and 3.5 //m. The complete photometric set of values 
with the observed and dereddened fluxes for each source is pre- 
sented in Appendix A. 

Most HAeBe stars are well-known variable stars, a problem 
that is relevant to our case because we used photometry recov- 
ered from the literature in the past decades. The variability af- 
fects the B-V color and produces some additional uncertainty in 
the extinction estimates. 



4. Model 

As mentioned in Sect. |3.U a two-component model is required to 
fit the observed SEDs. The disks were fitted usin g a passive, irra- 
diated circumstellar disk model described by lAlonso-Albi et all 
(2008) allowing us to consider different grain populations for 
the midplane and surface layers. Grain populations were char- 
acterized by a silicate/graphite mixture, the maximum grain size 
(o-max), and the slope of the grain size distribution (/?). The value 
of p was fixed to 3.5 and a standard grain mixture (SM: 86% 
silicate/14% graphite) were assumed for the midplane and sur- 



MWC 137 
R Mon 
MWC 1080 
MWC 297 
LKHa 215 
ZCMa 



SM a„ M _ t : 
SM a„ M . v : 
SM a„ M _ t : 
SM a„„„: 
SM a mflJ .: 
SM a,„„ v : 



= 1 cm 
= 1 cm 
= 1 cm 
= 1 cm 
= 1 cm 
= 1 cm 



SM a,„ 



= 1 yum 



86% graphite 



SM a,„ 
SM a„, 
SM a,„ 



100', 

50% 



= 100 yum 
= 1 yum 
= 100 yum 



60% graphite. 



100V 

50% 



= 1 yum 



= 1 yum 



The index of the power-law grain-size distribution is fixed (p=3.5). As 
a first step, we assumed the standard mixture (SM: 86% silicates/14% 
graphite) in the midplane and surface layers in all the sources. In the 
cases of R Mon and Z CMa, we needed to change the dust composition 
in the surface to fit the SED. The range of allowed values are indicated 
(small figures) in the Table. See Sect. 5.7 for a detailed discussion of 
the uncertainties. 



face layers. In two cases, R Mon and Z CMa, we were unable 
to obtain a reliable fitting with these assumptions and we needed 
to vary the grain mixture in the surface layer. These cases are 
discussed in detail in Sect. 5.7. 

In the fitting process, we consider that the emission at wave- 
lengths shorter than 10 /mi and the fluxes measured by interfero- 
metric techniques (PdBI, VLA) arise only in the star/disk system. 
This is based on the folowing scientific arguments: Firstly, the 
dust emitting at <10 /mi has a temperature of >300K. Assuming 
that the star is the only heating source in the region, the dust 
should be located at radii < 100-300 AU, which is the typical 
size of a circumstellar disk. Secondly, the spatial resolution of 
our interferometric mm and cm observations is usually below 
1000 AU. It is, again, reasonable to consider that most of the 
emission detected with the interferometers arise in the disk, al- 
though a contribution from the envelope may also exist in some 
cases. The emission observed at mid-IR and FIR wavelengths 
(Spitzer, SCUBA) are, however, assumed to be the sum of the 
disk and the envelope components. This is consistent with the 
MIPS images at 24 /mi, which show nebulous emission around 
MWC 1080, MWC 297, R Mon, and LkHa 215. 

The envelope was modeled using a ray-tracing code. The 
dust mass, density law, and geometry were determined by the 
fitting procedure. Our aim was to fit the observed emission and 
the visual extinction towards the star. Since the visual extinction 
was determined by the amount of mass along the line-of-sight, it 
constrains the geometry of the envelope. As discussed in Sect. [5] 
the morphology of the continuum submillimeter and millimeter 
maps were also taken into account in constraining the envelope 
geometry. We modeled each envelope by considering two possi- 
ble geometries, a sphere and a toroid. The dust temperature was 
assumed to be described by a power law, T(r) = T (r/r o y 05 , 
where the values of T and r depends on whether the enve- 
lope was shadowed or not by the inner circumstellar disk. In 
the non-shadowed case, typical of a sphere geometry, the effec- 
tive temperature and radius of the star were used as values for 
T„ and r . In the shadowed case, typical of toroidal envelopes, 
we extrapolated the dust temperature in the disk midplane. The 
du st opacities in the envelope we re assumed to be those reported 
bv lOssenkopf & Hennind(ll994l) for protostellar cores, using the 
standard MRN dust model with icy mantles. The resulting emis- 
sion was convolved by the beam of the observations before being 
compared with the observed SED. 

The goodness of fit was estimated us- 
ing a dispersion parameter (x) defined to be 
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Table 3. Summary of observational results 



Name 


3.6 cm 


1.3 cm 


7 mm 


2.7 mm 


1.3 mm 


Spec, index 1 


FpdBI 1 F 2 j C MT 




(mJy/beam) 


(mJy/beam) 


(mJy/beam) 


(mJy/beam) 


(mJy/beam) 


(1.3-2.6 mm) 


(1.3 mm) 


MWC 137 




0.9 


1.3 


4.1 


7.1 


0.79/- 


0.075 


RMon 




0.9 


1.3 


4.1 


11.8 


1.53/2.14 


0.15 


MWC 1080 








<1.7 


3.1 


>0.87/>1.02 


0.013 


MWC 297 




23 


29 


149 


175 


1.01 3 /1.31 3 


0.39 


LkHtt 215 








<0.6 


<1.5 


-/- 


<0.026 


ZCMa 


2.2 


2.1 


2.1 


8.5 


26 


1.61/1.88 


0.037 



1 The 1.3mm-2.7mm spectral index before (left) and after (right) subtracting the free-free emission. 

2 See references in Appendix A. 

3 Beam dependent, see Section [53l 



Table 4. Disk modeling results 1 



Name 


Dust 


Inner radius 


RimT 


Rim height 


Density 


Inclination 


Outer radius 


Dispersion 




Mass (M Q ) 


(AU) 


(K) 


(AU) 


slope 


(degrees) 


(AU) 





MWC 137 
RMon 
MWC 1080 
MWC 297 

LKHo-215 
ZCMa 



1.0x10- 



1-5 2xW- 4 
5.vlO- 6 

1.4x10---:: 

5.0xl0- 5 ^ 10 4 _ 5 
5.3 xlO- 4 

4.5*10 4 

4 0xl0- 4 60110 

-8 9xlQ- s 
4j10" 8 
4 1.4xl0" 3 
3.5.X-10- 4 



6.0x10 
7.0x10 



8.4 
18.2 
5.6 
5.0 

200 100 

iw 300 

3.5 
7.5 



1500 
1800 



1500 



1.0 
1.2 



1.9 



-0.45 °f. 
-0.55 :°| 

-0 77 -° 6 
u. / ; _ 9 



L - H -1.8 



-0.7 



>-1.0 



80?" 
60™ 
83?° 

5 ° 



30 f 
30 g 



18 25 
15 200 

77 100 
' ' 60 

28.5 30 
300™ 
9.2 20 

180 32 



0.20 
0.09 
0.08 
0.19 

0.1 
0.15 



1 In the case of fitted parameters, we provide in small figures the range of allowed values. These values correspond to numerical errors. See 
Sect. 5.7 for further discussion of the uncertainties in the disk parameters. 



Table 6. Envelope modeling results. 



Name 


Dust mass 


Inner radius 


Outer radius 


Density 


Inner temperature 


Geometry 


Inclination 




(Me) 


(AU) 


(AU) 


slope 


(K) 




(degrees) 


MWC 137 


0.005 


5000.0 


11700.0 


-1.0 


61 


Sphere 




RMon 


0.008 


700.0 


12000.0 


-1.4 


62 


Toroid 1 


70.0 


MWC 1080 


0.025 


6000.0 


12000.0 


-1.0 


44 


Toroid 


80.0 


ZCMa 


0.05 


2000.0 


5000.0 


-0.6 


42 


Toroid 1 


30.0 


LKHo-215 


0.0015 


3500.0 


7200.0 


0.0 


57 


Sphere 





1 The dust temperature in the envelope is assumed to be described by a power law, T(r) = Tj„(r/rj„) 03 . See Sections 4 and 5 for a more detailed 
discussion. 



X = VZ ((Fobserved ~ F modeled) I F observed) 1 1 n, where F is 

the flux and the sum is extended for all detections in the SED 
dominated by the disk emission from 2 /mi and beyond, and 
n the number of detections. In Table |4] we show the best-fit 
(minimum value of x) solution for each source. 

5. Individual sources 

5.1. R Mon 

The SED towards R Mon is the only in our sample for which 
the first peak occurs at 4-5 /mi instead of at 2 /mi (see Fig. [TJ. 
As a result, the SED can be fitted by a Chiang & Goldreich- 
like model without the presence of an inner rim. However, an 
inner rim cannot be discarded because we lack of photometry 
between 2 /mi and 3 /urn, and extinction at these wavelengths 
is very high. The near and mid-IR part of the SED is well fit- 
ted by the emission of a circumstellar disk with r OIi ,~150 AU. 
The interferometric 1.3 mm and 2.7 mm data allow us to con- 
strain the mass of dust and estimate the grain sizes. The dust 



mass in the disk around R Mon is found to be 1.4x10 4 M Q 
and grain growth has proceeded to s izes of ~ 1 c m in the mid- 
plane (see Table H and iFuente ft all d2003ll2006l) ). |Fuente et al.l 
d2006l) and lAlonso-Albi et al.l d2007l) demostrated the existence 
of a gaseous, flat disk around R Mon. The existence of large 
grains in the disk midplane is consistent with this geometry. 
Large grains produce a decrease in the dust opacity and a more 
rapid photoevaporation of the external layers of the gaseous disk. 

The difference between the single-dish and interferomet- 
ric fluxes at 1.3mm is clearl y ind icative of the presence o f an 
envelope. IClose et al] d 19971) and iMovsessian et all d2002l) in- 
ferred that the nebula illuminated by R Mon has an hourglass 
shape, and an inclination angle of ~20° with respect to the 
sky plane. The same kind of hourglass morphology i s observed 
in the 12 CO 1— »0 and 2— >1 images reported by IFuente et al.l 
(2006). To model the envelope, we used a toroid with the same 
inclination as observed and shadowed by the inner disk. The 
available photometry includes Spitzer mid-IR observ ations with 
a beam size of 30" and JCMT observations from Mannings 
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SED fitting for R Hon 




1 ID 1 10 2 10 3 10* 10 5 

Wavelength (/xm) 



Fig. 1. Left: Observed SED and model predictions for R Mon. Right: Sketch of our model. The different emission components that 
contribute to the SED as predicted by our model are drawn in different colors. The emission of the star itself is drawn in green. 
The disk emission is separated into three components: inner rim emission (magenta), surface layer emission (red), and the midplane 
emission (blue). Note that the disk around R Mon is fitted without the presence of an inner rim. The emission from the envelope 
is drawn in yellow. The free-free emission is shown in magenta. The disk emission and the disk+envelope emission appear as two 
continuous, black lines. 



(1994) with a beam size of around 18". We extended the enve- 
lope to cover the Spitzer beam size, using a temperature profile 
T(r) = 62(r(At/)/700r a5 K (see Table|<§. The mass of dust in 
the envelope was found to be 0.008 M . The density profile of 
the envelope was described by a power law with slope -1.4. This 
means that an important fraction of the mass is close to the inner 
radius of the envelope, which is found to be 700 AU. The dust 
temperature at this radius is 60 K. 



5.2. Z CMa 

Z CMa is a binary system with a HBe star and a FU Orionis 
companion. The most recent estimate of the stellar masses of 
Z CM a and its companion were 16 M and 3 M , respec- 
tively (Ivan den Ancker et alj 120041) . The available photometry 
only provides the combined emissi on of both componen ts, with 
the exception of the observations by Koresk oet al.l (Il99ll) which 
resolved the two stars. The optical observations shown in Fig. 
|2] should be considered to be a combination of both com- 
ponents. We fitted the optical observations by assuming that 
the HBe star contributes ap proximately 30% of the total flux 
(I van den Ancker et al.ll2004l) . 

The SED of Z CMa is shown in Fig [2] The shape from near 
to far IR is flat with an abrupt slope in the emission beyond 100- 
200 fim. These two characteristics indicate that this is a massive, 
compact object. Obviously, the disk properties are difficult to 
derive and the existence of two circumstellar disks, one around 
each star, and/or a circumbinary disk cannot be discarded. We 
attempted to reproduce the NIR, mid-IR and FIR SED by assum- 
ing one single disk around the HBe star. We were unable to fit the 
IR SED by assuming a value of 16 M for the stellar mass, but 
we obtained a reasonably good fit by assuming a stellar mass of 
12 M which is still consistent with being a Be star. On the basis 
of mainly our interferometric mm observations, we derive a mass 
of ~7.0xl0~ 4 M for the circumstellar dust. The 1.3mm/2.7mm 



spectral index was fitted most accurately with grains of maxi- 
mum size ~ 1 cm. 

The envelope around Z CMa is the most massive and com- 
pact in our sample (see Fig. [2] and Table [6j, which may be in- 
terpreted as proof of the youth of this object. It is well fitted by 
a shadowed toroid of r,„= 2000 AU, r ouf =5000 AU, and an in- 
clination angle of 30°. The abrupt slope in the sub-mm region 
implies that the envelope cannot be extended beyond an outer 
radius of 5000 AU. For this source, some measurements derived 
for different b eams of fluxes at 450 yum, 850 yum, and 1.3 mm 
are available dDent et alJll998t ISandell & Weintraubll200Tl) . By 
fitting the observed flux by different beam models, we inferred a 
density slope of -0.6 for the envelope, a value close to the density 
slope of -0.7 found in the disk. 

The derived extinction from the envelope model (Table |A~8l ) 
is since the inclination of the toroid is small, 30°, close to a 
pole-on view, and the star itself is visible through the conical 
cavity excavated by the outflow. However, by increasing the in- 
clination to 50°, the extinction increases to a value a factor of 
10 higher than the observed value. An outer cold, spherical en- 
velope could possibly account more accurately for the observed 
extinction. 

5.3. MWC 297 

At a distance of only 250 pc, MWC 297 is the closest star in our 
sample. It is highly reddened with a B-V color of 2.0, with no 
known binary companion. At optical wavelengths, a value of R v 
= 3.5 accounted for the extinction better than the value of R v = 
3.1. At 2 //m, the SED is almost flat with a small prominence 
that could be interpreted as an inner rim (see Fig. [3]). However, 
this feature is unclear. We obtained reasonably good fits by us- 
ing a face-on disk (i ~ 5°) model with an inner rim, or a highly 
inclined disk (i ~ 80°) witho ut an inner rim. O n the basis of IR 
interferometric observations, Ack e et aTl d2008l) proposed the ex- 
istence of a low inclination disk (z less than 40°) with dust closer 
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Fig. 2. Up: Z CMa maps obtained with PdBI. Contours in the 1 .3mm image are represented in steps of 0.77 mJy/beam, starting from 
2.31 mJy/beam. In the 2.6mm image, the contours are from 1.8 mJy/beam, in steps of 0.6 mJy/beam. Down-Left: Z CMa map at 
1.3 mm obtained with the 30m bolometer. Contours are represented in steps of 20 mJy/beam starting from 50 mJy/beam (3<x level). 
Maximum intensity is 293 mJy. The rectangle indicates the 20 arcseconds side region shown in the upper panels. Down-Right: 
The observed SED and our model predictions for Z CMa. The different emission components are drawn following the color code 
described in Fig.Q] 



to the star than the dust sublimation radius. We selected the low- 
inclination solution to be con sistent with pre vious near-IR inter- 
ferometric measurements bv lAcke et al.l d2008l) . Since our model 
is based on hydrostatic equilibrium, it cannot account for the ex- 
istence of dust closer to the star than the sublimation radius. In 
our model, we allow the rim temperature (Trim) to vary between 
1500 K (silicates sublimation temperature) and 2000 K (graphite 
sublimation temperature) and the inner radius is determined by 
T rim . The fact that the flux decreases at wavelengths longer than 
10 /im implies that the circumstellar disk should be small. The 
most consistent fit is obtained with T nm ~1800 K, an outer radius 
of only ~28 AU, a density slope of ~-0.77, and an inclination 
angle of 5° + 5° . This small disk is optically thick at 1.3mm 
and, consequently, the dust mass and grain properties are not de- 
termined (see Table|5]l. We assume a max = 1 cm in the midplane. 

The slope of the SED at mm and sub-mm wavelengths is 
inconsistent with that expected from an envelope with standard 
interstellar-medium (ISM) grains. The 1.3 mm/2.7 mm spectral 



index is also too low to be adjusted by thermal dust emission. 
The problem could be the different beam sizes of the interfer- 
ometric observations at 1.3 mm and 2.7 mm. In the case that 
the emission at 1.3 mm and 2.7 mm is more extended than the 
beams, the measured 1.3 mm/2.7 mm spectral index would be af- 
fected by the different angular resolutions. To confirm this pos- 
sibility, we compared our 1.3 mm flux with the previous mea- 
surement carried o ut with the SMA providing a beam of ~3" 
(iManoi et al.l2007h . Our flux 58% of the flux values measured by 
iManoi eta l. (2007), indicating that our 1.3mm observations re- 
sol ve the 1 ,3mm emission. Thus, we adopted the flux measured 
by IManoi et al.J d2007l) to derive the 1.3mm/2.7mm spectral in- 
dex. Even using this value, the derived spectral index was too 
low to be due to the ISM grains expected in the envelope. The 
spectral index was, however, more similar to that found in cir- 
cumstellar disks, where grain growth occurred. This prompted 
us to consider the existence of a second disk component to- 
wards MWC 297. This second disk could be due to a binary 
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companion, or a ring around the main disk. Since no binary 
companion was detected close to this star (at a distance of ~ 
0.7"or 180 AU), the second option is more likely. We found 
good agreement with observations by assuming a ring of large 
grains (~ 1cm), a dust mass of 4xl0~ 4 M Q , an inner and outer 
radius of 200 AU and 300 AU respectively, and a temperature 
profile T(r) = 125(r(At/)/200r 05 K. 

The JCMT map (bottom-left panel in Fig. [3]) indicates that, 
in addition to the small disk and ring, there is a more extended 
envelope surrounding this object. Nebulous emission is also ob- 
served in the MIPS 24 fim image, suggesting that the star is still 
embedded in the molecular cloud. However, most of the emis- 
sion at millimeter wavelengths (92%) indicates a small region 
(~3") around the star (see Appendix A. 7). This is consistent with 
our two-disk-component fit. 

5.4. MWC 1080 

MWC 1080 was classifi ed as a BO star. Its distance e stimates 
ranged f rom 1000 pc (lEisner et al.ll2004t iFuente et alll2003l) to 
2100 pc (lAcke et al.ll2005l) . We adopted the distance of 1000 pc, 
since it is far more consistent with the observed stellar emission 
after the extinction correction. 

The SED towards MWC 1080 has a clear bump at 2 yum, in- 
dicating the presence of a rim with T f /~1500 K (see Fig.[4]l. The 
dust temperature in the inner rim was determined from the peak 
at 2 /urn. The mid-IR part of the SED was well fitted with a disk 
of r oi „~77 AU and an inclination angle of ~80°. This inclination 
angle is totally consistent with MWC 1080 being an eclipsing 
binary star, for which the amplitude in the var iation of the ap- 
parent magnitude is wide (Gran kln et alj|1992l) . The dust mass 
that better accounts for the 1.3 mm interferometric flux is found 
to be 5xl0~ 5 M . The maximum grain size is not determined 
since we have not detected the disk at 2.7mm. We obtain a good 
fit for grains of ~1 cm in the midplane. 

The existence of an envelope is certain since the observed 
1.3mm interferometric flux is a small fraction (1.3%) of the to- 
tal flux measured with the JCMT. The interferometric 1.3 and 
2.7mm images also reveal the existence of several clumps sur- 
rounding the circumstellar disk (see Fig. [4] and Fuente et al. 
2003). On the basis of the morphology observed in our interfero- 
metric images, we used a toroid with an inner radius of 6000 AU 
(6" at the distance of MWC 1080) to fit the envelope. The best 
fit solution for the envelope is a non-shadowed toroid of an outer 
radius of 12 000 AU (12"). Our model fits all the SED points ex- 
cept the IRAS flux at 100 fim. This is unsurprising since the wide 
IRAS beam at 100 fim (~2 r ) probably detects emission from the 
foreground molecular cloud. 

The extinction estimate from the B-V color index is 40% 
higher than that predicted by our envelope model. This is also 
consistent with the scheme of the star+disk+toroid system im- 
mersed in a lower-density molecular cloud. 

5.5. MWC 137 

The SED of MWC 137 does not exhibit any evidence of disk 
emission at mm wavelengths. All the mm-cm fluxes can be fit- 
ted by a single component of spectral index a = +0.76+0.01. 
This spectral index is consistent with that expected in the free- 
free emission arising in the stellar wind. Although a value of 
a =+0.6 is expected for an ionized isotropic wind, small devi- 
ations of this value can be explained by a different geometry 
or a partially ionized wind. However, some excess is visible at 
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Fig. 5. The observed SED and our model predictions for 
MWC 137. 



SED fitting for LKHa 215 




Wavelength (jum) 

Fig. 6. The observed SED and our model predictions for 
LKHa 215. 



IR w avelengths, which i s suggestive of the presence of a small 
disk. lFuente et al.l d2003) derived an upper limit to the dust mass 
of 7x1 0~ 5 M assuming that the disk is optically thin at 1.3mm, 
the free-free emission spectral index is +0.6, and standard val- 
ues for the dust opacity and mean dust temperature (see Sect. 6). 
However, this value is quite uncertain since the disk is probably 
optically thick and the grain emissivity is unknown. We fitted the 
NIR part of the SED and obtained the best-fit solution (consis- 
tent with our mm observations) for a small disk of r ou ,~18 AU 
and a dust mass of 10~ 5 M . The absence of a large disk ensure 
that the envelope is prominent in the SED. The envelope is mod- 
eled as an sphere extended between 5000 and 1 1 700 AU, with a 
dust mass of 0.005 M (see Fig. [5j>. 
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Fig. 3. Up: MWC 297 maps obtained with PdBI. Contours in the 1.3 mm image are represented in steps of 4 mJy/beam, starting 
from 12 mJy/beam. In the 2.6 mm image the contours are from 6.6 mJy/beam, in steps of 2.2 mJy/beam. Down-Left: MWC 297 
map at 850 fim obtained with JCMT. Contours are represented in steps of 75 mJy/beam, starting from 75 mJy/beam. The small 
rectangle marks the 6 arcseconds side region shown in the upper panels. Down-Right: The observed SED and our model predictions 
for MWC 297. Note that in this case there is no envelope but two disk components. 



5.6. LKHa 215 

We did not detected this source at mm and cm wavelengths. This 
implies an upper limit to the dust mass of 9.0xl0~ 5 M , assum- 
ing optically thin emission and standard values of the dust tem- 
perature and emissivity. We modeled the entire SED to improve 
the upper limit to the mass of dust and to determine the envelope 
component (see Fig. [6). Our best fit, still although uncertain, so- 
lution is obtained for a disk with an outer radius of ~10 AU 
and a mass ~6xl(T 8 M . The envelope is described well by a 
dust mass of about -0.0015 M , located between ~3500 AU and 
-7200 AU. 



5. 7. Uncertainties in the disk mass and maximum grain size 

The dust mass and the grain properties determined by fitting the 
SED are affected by important uncertainties. The most impor- 
tant assumption is that the dust emission is optically thin at mm 
wavelengths. This is the case for typical circumstellar disks of 



rout ~ 100 AU and dust mass < 10 3 M Q . Only in the case of 
small disks, the 1.3mm emission could become optically thick. 

We discuss R Mon as a representative case of an optically 
thin disk. Our best fit solution is for a disk with r o „,= 150 AU 
and a dust mass of 1.4xl0~ 4 M . This disk is optically thin 
at 1.3mm. A disk of this size needs to have a dust mass of 
~0.01 M , i.e. a factor of 100 higher than our estimate, to be- 
come optically thick (tj 3„, m ~1). But in this case, the flux at 
1.3mm would be more than one order of magnitude higher than 
the measured flux (see Fig. [7]). In the optically thick case, the 
flux at 1.3mm depends only on the area of the disk projected on 
the plane of the sky and the dust temperature across the disk sur- 
face, Fi 3„ lm ~ BJT ' s )7irl ut cos(i). The only way of reconciling 
the predicted mm flux with observations is by reducing the size 
of the disk and/or changing the disk inclination. But in that case 
we obtain a deficient emission in the NIR and mid-IR range as it 
is illustrated in Fig. [7] The same argument is valid for Z CMa. 

A different case is the small disk surrounding MWC 297. In 
this case, the disk is optically thick and the derived mass is only 
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Fig. 4. Up: MWC 1080 maps at 1.3 mm and 2.6 mm obtained with PdBI. Contours in the 1.3 mm image are represented in 
steps of 0.45 mJy/beam, starting from 1.35 mJy/beam. In the 2.6 mm image, the contours are from 0.45 mJy/beam, in steps of 
0.15 mJy/beam. Down-Left: MWC 1080 map at 1 mm obtained with the 30m bolometer. Contours are represented in steps of 
10 mJy/beam starting from 40 mJy/beam (3cr level). Maximum intensity is 95 mJy. The 26 arcseconds side square indicates the 
region shown in the PdBI 2.6mm map. Down-Right: The observed SED and our model predictions for MWC 1080. 



a lower limit to the true disk mass. However, it is unreasonable 
to propose that the dust mass in the small and, probably, evolved 
disk around MWC 297 is larger than that within the less evolved 
disk around R Mon. One would expect that the mass of dust 
and gas in the disk decreases because of photoevaporation (see 
Sect. 6). Thus, our mass should be accurate to within an order of 
magnitude. 

The SED towards MWC 1080 is well described by a partially 
optically thin disk as shown in Table|4] But it can also be reason- 
ably well reproduced by an optically thick disk with r OM( = 50 AU 
and a dust mass > 10~ 3 M . Following the same argument as for 
MWC 297, it is unreasonable to assume that the dust mass in 
this small disk is one order of magnitude larger than in the less 
evolved disk around R Mon and similar to those found in TTs. 
Thus, we keep the optically thin solution, although we are aware 
of the uncertainty in our analysis. 

In the case of optically thin disks, we can learn about the 
grain properties from the 1.3mm/2.7mm spectral index. In the 
Rayleigh-Jeans region of the spectrum, the dust emission is oc 



y (2+/s) wn ere is the opacity index that can be accurately de- 
rived by fitting the observed SED. The value of /3 in the submm- 
mm ran ge is a n excellent indicator of the grain size distribution 



(Draine 2006). The grain size distribution is usually taken to be, 
n(a) = noa~ p for a < a max , where no is a normalization factor, 
p is a parameter in the range 2.5-3.5, and a max is the maximum 
grain size. The value of f3 is between 1 and 2 for small grains 
(sub-micron-sized grains) but becomes lower than 1 when a max 
is a few mm. To have a value of j6<0.5, we need grains with 
a m ax~ 1 cm. This is the case for R Mon and Z CMa. However, 
there is some degeneracy between the index of the grain size 
distribution, p, and the value of the maximum grain size, a max . 
A slope p=2.5 means that a larger number of grains are closer 
to the maximum grain size than for a more typical slope of 
3.5. By assuming p=2.5 we would reach the same values of j3 
with grains slightly smaller. Even in this case, we would require 
grains with radii of nearly 1 cm, i.e. significant grain growth 
must have occurred in the disk midplane to have the values of 
<0.5 that we obtain in our sources. We cannot determine grain 
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Fig. 7. Left: The SED of R Mon cannot be fitted assuming an optically thick disk. In this case, we obtain a large excess at mm 
wavelengths (green short-long dashed line). To reduce this excess, we can adjust the other disk parameters, such as the inclination 
angle or the outer radius, to reduce the effective surface of the disk. But in this case, we obtain a deficient emission at near and 
mid-IR wavelengths (blue dashed line). The red continuous line shows our optically thin solution. The lower-left chart shows the 
evolution in the dispersion parameter (x). Right: The SED of R Mon and Z CMa cannot be fitted using a standard grain mixture (86% 
of astronomical silicate and 14% of graphite). The pink line shows our best fit for a non-standard mixture with 86% of graphite. 
With 100% of graphite the fit is still excellent, but slightly worse. When the fraction of silicates becomes important (50% or higher) 
the model is inconsistent with the observations. 



sizes larger than a few cm on the basis of mm observations. The 
1.3mm/2.7mm grain opacity index tends asymptotically to for 
O-max > a f ew cm - in the case of optically thick disks, we cannot 
determine the grain size from the SED. We assumed a max ~ \ cm 
in the midplane for MWC 1080, MWC 297, MWC 137, and 
LkHa 215, since grain growth is expected to have proceeded in 
these more evolved disks. 

The dust composition in the midplane cannot be inferred 
from the SED. For this reason, we assumed the standard mix- 
ture in all of our sources (see Table [5]). The silicate feature at 
~9.8 //m provides, however, some information about the dust 
composition in the surface layer. In all of our sources, the sili- 
cate feature is weak or absent. The absence of this feature does 
not imply the absence of silicate grains but the lack of silicate 
grains at a temperature of ~800 K. The silicate grains could be 
either too far from the central source and thus not heated to tem- 
peratures sufficiently high to allow emission in the mid-IR, or 
these grains could be too large to be heated efficiently. To inves- 
tigate the second possibility, we tried to fit SEDs by assuming a 
standard dust composition and varying the a max in the surface be- 
tween 0. lyi/m up to 100 fim. We were able to reproduce the disks 
around MWC 137, MWC 1080, MWC 297, and LkHa 215 (see 
Table [5]), but failed to account for R Mon and Z CMa. In these 
cases, we needed to vary the dust composition to fit the mid-IR 
part of the SED. The most extreme case was R Mon in which 
the best fit solution was found for 86% of graphite in the surface 
layer (see Fig. [7]). This result is, of course, limited by the sim- 
plicity of our model. A different disk geometry is expected in 
a thermal equilibrium disk. Alternatively, the dust composition 
and/or grain growth could vary as a function of the distance from 
the star. Determining the dust composition of the disk surface is, 
however, beyond the scope of this Paper, and, in addition, does 



not affect our results because the disk mass depends only on the 
dust composition and grain size in the midplane. Therefore, we 
do not discuss this point further. 

The dust composition and maximum grain size in the 
midplane determine the value of dust emissivity at 1.3mm, 
Ki3mm, and consequently has an influence on the derived dust 
mass. Assuming a standard dust mixture and a max -\ cm, 
Ki.3mm=0.88 g -1 cm 2 . This value is close to the canonical value 
of 1 g~' cm 2 . For reasonable values of a mfl v in disks, i.e. of be- 
tween 1 mm and 1 cm, the value of Ki.3 mm varies by a factor of 
2-3 from the canonical value. This is the uncertainty in the disk 
mass estimate due to the uncertainty in the grain size in the mid- 
plane. 

Summarizing we conclude that our disk-mass estimates are 
accurate to within a factor of 2-3 in the optically thin disks 
around R Mon and Z CMa, and to within an order of magnitude 
for the optically thick disks around MWC 1080 and MWC 297. 

NIR interferom etric observations bv lEisner et alj d2007h and 
lAcke et alj d2008l) revealed the existence of hot dust in a region 
at < 1 AU from the star towards MWC 297 and MWC 1080. We 
are aware that our model is too simple to determine the prop- 
erties of the disk in a few AU around the star. First of all, in 
our model, the inner radius of the disk is determined by the dust 
sublimation temperature that we assume to be between 1500 K 
(silicate sublimation temperature) and 2000 K (graphite subli- 
mation temperature). We do not consider the possib ility of the 
existen ce of dust closer to th e star, as proposed by lEisner et all 
d2007l) and lAcke etail d2008l) . Secondly, this extremely hot dust 
(at less than 1 AU) would emit radiation at <2 /mi and the SED 
at these wavelengths is dominated by the emission of the stellar 
photosphere. Our model and method (SED fitting) are indeed in- 
adequate to determine the properties of the innermost regions of 
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the disk (r<a few AU). NIR interferometric observations capa- 
ble of separating the star from the disk are required to probe this 
region. 

We conclude that the SED fitting technique is sufficiently 
robust to reproduce the global parameters (dust mass, grain 
growth) of the outer disks. We emphasize that our analysis de- 
scribes reasonably well the SED for a wide range of wavelengths 
(3 jt/m-3.6 cm) with a consistent set of star, disk, and envelope 
parameters. Our results are also consistent with the morphology 
observed in the millimeter and mid-infrared images. 



6. Discussion 

Our cm and mm interferometric observations have provided 
valuable insight into the structure of the circumstellar matter 
around HBe stars. The first result is that HBe stars are usually 
surrounded by a massive envelope. The mass of the envelope 
is always larger than the mass of the disk and dominates the 
emission at FIR and millimeter wavelengths. Extended emission 
arising in the envelope is also observed in the MIPS images at 
24 //m, which shows that the envelope emission significantly 
contributes to the SED even at mid-IR wavelengths. The exis- 
tence of a nebula w as also inferred for Spit zer images towards 
the AO star VV Ser dPontoppidan et al.ll2007h . Our observations 
corroborate that interferometric observations are required to de- 
rive the properties of the disks around these massive stars, and a 
disk+envelope model is essential when interpreting the SED. 

We have different types of disks in our sample. Towards 
Z CMa and R Mon, we detected disks with sizes of r ou/ >100 AU 
and dust masses of 10~ 4 -10" 3 M . These values do not differ 
from the estimated sizes of disks around TTs and HAe stars. 
Both of these disks were detected in the CO J=l— >0 and 3-2— >1 
rotational lines (Fuente et al. 2006, Alonso-Albi et al. 2009). 
Z CMa is a controver sial case. Its accretion rate of 10~ 3 M /yr 
dHartmann et alj 19891) is one of the highest known around a pre- 
main sequence object, and Z CMa is surrounded by a compact, 
massi ve envelope and powers an energetic outflow dPoetzel et alj 
119891) . It is probably the youngest object in our sample, and 
it probably contains the most massive disk. However, the dust 
emission could originate in a circumbinary ring and/or a disk 
associated with the FU Orionis companion. 

The circumstellar disk around MWC 1080 is smaller 
(r o »f~80 AU) and less massive (5xl0~ 5 M ) than those around 
Z CMa and R Mon. However, the envelope contains a signifi- 
cant amount of dust, which suggests that the star is still quite 
young. This disk has not been detected in the CO J=l — >0 and 
2— >1 lines (Alonso-Albi et al. unpublished data). 

MWC 297 is a special case because its disk consists of 
two components. The inner disk has a radius of 28 AU, simi- 
lar to the more evolved sources discussed below, while the outer 
ring has a radius of 200 AU, more similar to the young disks 
around R Mon and Z CMa. The disk has not been detected by its 
CO J=l — >0, J=2— >1 (Al onso-Albi et al. unpublished data) nor 
J=3->2 dAcke et al.ll2008l) rotational line radiation. 

Although the mass of dust in the circumstellar disk 
varies between one source and another (10 5 — 10 3 M ), the 
1.3mm/2.7mm spectral index is low (<2.2) in all sources. This 
1.3mm/2.7mm spectral index indicates that grain growth has 
proceeded to sizes of 1 cm even in the youngest sources of 
our sample, R Mon and Z CMa. As discussed later, early grain 
growth has important consequences for the evolution of the 
dusty disk. The non-detection at mm wavelengths of the disks 
around MWC 137 and LkHa' 215 also implies the existence of 



compact (<20 AU) and/or presumably light disks (<0.001 M ), 
if any, around these stars. 

6. 1 . Disk masses around Herbig Ae/Be stars 

To investigate the possible dependence of the disk mass on the 
stellar mass and age, we compiled a list of disks detected at 
mm wavelengths (1.1-2. 7mm) around HAeBe stars (see Table 
IA. 13b . It is important to calculate the disk masses in a uniform 
way. We calculated all disk masses, including for our own sam- 
ple, by assuming optically thin emission at mm wavelengths, 
an av erage disk tempe rature that depends on the stellar spectral 
type (Nat ta et al.l2000l) . a dust opacity Ki.3 mm =1.0 cm 2 g -1 , and a 
dust spectral index [} — 1.0. We subtracted the free- free contribu- 
tion in the sources previously studied by our team dFuente et al.1 
I2003L 120061: lAlonso-Albi et alj 120081) . Then we calculated the 
disk masses by assuming that all the 1.3mm/2.7mm emission 
originates in the dusty disk. To estimate the uncertainty in the 
disk masses derived in this "simple" way, we compared the disk 
masses for our sample with those previously derived by model- 
ing the full SED in this Paper. The values agree within a factor 
2-3, which suggests that our method is robust to within the sta- 
tistical uncertainties. 

We also calculated the ma ss and age of all o f these stars us- 
ing the evolutionary tracks of ISiess et al.l (HOOO) and the values 
of luminosities and effective temperatures found in the literature 
(see Table lA~T3b . In some cases, the luminosities and effective 
temperatures were mod ified slight l y (w ithin the expected uncer- 
tainty) to fit the closest lSiess et alj d2000h evolutionary track and 
the spectral type of each star. Stars with masses higher than 7 M 
do not undergo a pre-main sequence phase and their ages cannot 
be estimated using the HR diagram. In these cases, we evalu- 
ated the age on the basis of different arguments. For Z CMa, 
we assumed an age of 0.05 Myr consistent with the time re- 
quired to form a star of 16 M with a mean accretion rate of 
3x10 4 M yr _1 . This star could have been even younger if the 
accretion rate had remained close to 10~ 3 M /yr throughtout 
the formation process. MWC 137, MWC 297, MWC 1080, and 
R Mon are very likely young main-sequence objects. For these 
objects, we assumed that the age was lower than 1 My r based on 
the post-main sequence tracks of Schaller et al. ( 1992). For clar- 
ity and to differentiate them from other sources in which we have 
been able to estimate stellar ages, we plot them in the 0. 1 Myr 
slot in the disk mass versus stellar age diagram. 

In Fig. [8] (right panel), we present the total disk masses 
(gas+dust) as a function of the stellar mass assuming a gas/dust 
ratio of 100. We can easily distinguish three regions in the plot. 
Firstly, for stars with masses of <3 M , there is a cloud of points 
around a mean value of ~0.04 M . However, there are a few 
stars with light disks (<0.01 M ). Secondly, in stars with masses 
of 3-7 M , there is a group of massive disks (>0. 1 M ). Most of 
these disks masses were derived from 2.7mm observations and 
were probably overestimated because of the unsubtracted free- 
free emission. In the case of HD 142527 and VY Mon, the disk 
masses were derived from single-dish observations and we can- 
not discard the contribution of an associated nebula. Thus, these 
massive disks requires further confirmation. Thirdly, in more 
massive stars (>7 M ), mainly our entire sample of HBe stars, 
the disk masses decrease by a factor 5-10 and become <0.01 M . 
The only exception is Z CMa, a particularly young object as dis- 
cussed above. Since it is a binary containing a 16 M and 3 M 
stars, the circumstellar disk could be associated with the lower 
mass component and/or a circumbinary disk. Another interest- 
ing case is MWC 297. If we consider only the inner disk, it is 
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Disk mass as function of star age Disk mass as function of star mass 




10 5 10 6 10 7 10' 

Star age (yr) Star mass (M ) 



Fig. 8. Compilation of disks masses (gas+dust) as function of stellar age {left) and the stellar mass (right). Our observations with 
black circles. In blue diamonds we represent the disk masses derived from 1.3 mm interferometric observations from other authors. 
With green inverted triangles we represent the disk masses derived from 2.7 mm interferometric observations. In these stars the disk 
mass is probably overestimated because of the possible (and not subtracted) contribution of the free-free emission. The red circles 
correspond to single-dish 1.3 mm observations. The disk mass in these cases could also be overestimated because of the existence 
of an envelope. 



located in the region of the "light disks" (<0.01 M Q ) related to 
massive stars. If we consider the sum of the two disk compo- 
nents, the disk mass is similar to those of the disks associated 
with HAe stars. 

In Fig. [8](left panel), we show the disk masses as a function 
of the stellar age. We can also distinguish the three regions de- 
scribed above. The cloud of points of stars with masses <3 M 
and disk masses ~0.04 M corresponds to stars of ages between 
a few Myr and 10 Myr. The most massive disks corresponds to 
stars of masses between 3-7 M Q and ages of around 1 Myr. The 
circumstellar disks around stars with masses of >7 M exhibit a 
large mass dispersion, ranging from >0.1 M to <0.01 M , i.e., 
two orders of magnitude, although mainly at the lower end of 
this range. 

To attempt to understand this behavior, we recall that we 
only consider pre-main sequence (visible) objects. Thus, there 
is a lower limit to the age of the star that corresponds to the 
birthline and depends on the stellar mass. The upper limit to 
the stellar ages for our disk detections is given by the disk life- 
time which probably also depends on stellar mass. For Herbig 
Ae stars with masses between 1-3 M , this period (from the 
birthline until the disk disperses) corresponds to stellar ages 
of between a few Myr and 10 Myr. During this phase the disk 
mass does not vary significantly with stellar age. For stars with 
masses in the range 3-7 M , this phase lasts from 10 5 yr to 
~2 Myr. For the most massive stars, >7 M Q , once the star is vis- 
ible the disk dispersal is so rapid that instead of having a cloud 
of points around a mean value, we have a large dispersion of 
the disk masses (2 orders of magnitude), and most of objects 
in the lower end (M<0.01 M ). This indicates that, although 
massiv e disks are found towards young, intermed iate-mass pro- 
tostars dCesaroni et al.ll2005t ISchrever et al.ll2006h and possibly 
Z CMa, after 1 Myr massive stars have already dispersed most of 



the circumstellar material and are surrounded by disks of masses 
lower than 0.01 M . 



Table 7. Dissipation timescale. 



Source 




Mass loss rate 


Life time 




(AU) 


(M /yr) 


(yr) 


MWC 137 


124.6 


8.9xl0~ 6 


1.1x10 s 


LKHff 215 


53.4 


3.4xl0~ s 


2.9x1 5 


RMon 


71.2 


4.1xl0- 6 


2.4xl0 5 


ZCMa 


106.8 


6.1xl0- 6 


1.6xl0 5 


MWC 297 


80 


6.0xl0- 6 


1.7xl0 5 


MWC 1080 


89 


4.6xl0~ 6 


2.2xl0 5 



1 Grav itational radius calculated following expression (1) by Ale xander! 
d2008h 

2 Time required to disperse a 1 M G disk. 



6.2. Disk photo-evaporation in massive stars 

Different dynamical processes contribute to the dispersal of gas 
in circumstellar disks: viscosity, winds and jets, dynamical en- 
counters with nearby stars, and photoevaporation by the central 
star and external illuminatio n. The relative importan ce of these 
processes was discussed by iHollenbach et all (120001) . who con- 
cluded that the dominant processes affecting disk evolution are 
viscosity in the inner regions (<10 AU) and photoevaporation at 
larger radii. Photoevaporation affects firstly the outer region of 
the disk, where the gravitational field is lower and the thermal 
heating can disrupt the disk (see e.g. the review by Alexander, 
2007). We can define a characteristic radius R g beyond which 
material evaporates and a gap opens in the disk. As soon as 
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the gap opens, the direct photoevaporation disperses the gaseous 
disk from the inside to the outside in a few 10 5 yr. 

While the gas evolution is outlined above, the dust evolution 
is mo re difficult to predict and depends on the maximum grain 
size (Garaud 2007). Particles are coupled to the gas by means of 
the frictional drag that depends on grain size. If the particles are 
far smaller than the mean free path of the gas molecules (A m f p ), 
gas and grains are coupled. Assuming typical gas densities, the 
grains with sizes smaller than 1 fim remain coupled to the gas. 
Thus, if grain growth is ignored, the evolution of the solid par- 
ticles relates to that of the gas and the dusty disk is rapidly dis- 
persed once the gap is opened. However, if grain growth is taken 
into account the evolution of the dusty disk is different. Although 
small particles are blown away with the gas, the large particles 
remai n forming a compact, dusty disk. lAlexander & Armitagd 
(120071) showed that there is a radius beyond which all solids 
are entrained in the gas. This radius is identified with the region 
where the maximum grain size is between 1-10 cm. 

We propose that the presence of compact, dusty disks around 
most HBe stars is the consequence of rapid disk photoevapora- 
tion in these massive stars. Since the ionizing flux in HBe stars 
is higher than in the cooler HAe and TT stars, the timescale for 
disk photoevaporation is shorter. In Table [7] we show the char- 
acteristic gap radius and the phot oevaporation rate for the stars 
in our sample calculated following Alexander (2008). The mass- 
loss rate is about an order of magnitude higher in Be stars than 
in Ae stars. While VV Ser as an AO star needs about 9 Myr to 
dissipate 1 M disk mass, HBe stars disperse their gaseous disk 
on a few 10 5 yr. This time is so short that the star could dissipate 
the outer part of the disk before becoming visible. As soon as the 
outer disk is dispersed, the entire gaseous disk is dispersed on a 
short timescale. This is consistent with the lack of CO detec- 
tion in the small disks, MWC 1080 and MWC 297, of our sam- 
ple. One important result of our mm-wave analysis is that grain 
growth begins early in disk evolution, and grains have grown to 
sizes of ~ 1 cm even before the star becomes visible. Thus, once 
the gaseous disk is dispersed, a compact, dusty disk formed by 
large grains remains. An interesting case is MWC 297 in which 
we find a small, circumstellar disk (r ol „=29 AU) and a remnant, 
outer ring with r,„=200 AU and r oi „=300 AU. This second ring 
could be the consequence of the decrease in gas density outwards 
from the star centre. The grains entrained in the retreating gas are 
left behind when the gas density is too low for the drag force to 
push them away. MWC 297 could be an object in-between the 
young disks associated with R Mon and Z CMa and the compact, 
dusty disk surrounding MWC 1080 and presumably, MWC 137 
andLkHo-215. 

Additional support for the rap i d ph otoevaporation of disks in 
HBe stars is found bv lBerne etall d2008l) on the basis of the PAH 
and small grains (VSG) emission in the surroundings of HBe 
stars. They suggested that the PAH/VSG emission from disks 
around B stars is absent. The observed emission emanates from 
the PDRs formed in the remnant of the parental cloud. They in- 
terpret this lack of emission as a sign that PAHs/VSGs have been 
photodissociated by hard UV photons in the disks. Since PAHs 
represent most of the UV optical depth, we conclude that this is 
also in good agreement with a scenario of photoevaporation in 
the disk around massive stars. 



7. Conclusions 

We briefly summarize the main conclusions of our study: 



- We have carried out a search for circumstellar disks 
around HBe stars using the NRAO Very Large Array 
(VLA) and IRAM Plateau de Bure (PdB) interferometers. 
Interferometric observations are necessary to separate the 
disk from the surrounding envelope, although the disks 
themselves usually remain unresolved. Our goal is to inves- 
tigate the properties of the circumstellar disks around inter- 
mediate mass stars to determine eventually their occurrence, 
lifetime, and evolution. Thus far, we have observed 6 objects 
with 4 successful detections. 

- Our first result is that the disk mass is usually only a small 
percentage (less than 10%) of the mass of the entire enve- 
lope. There are significant variations in the disk mass from 
source to source. Two disks of our sample, R Mon and 
Z CMa, have similar sizes to those found in TT and HAe 
stars. Z CMa is a FU Orionis object and cannot be directly 
compared with the others. The mass of the other large disk, 
R Mon, is a factor of 5 below that of an average HAe star. 
The disk around MWC 1080 is smaller 0w<100 AU) and 
less massive. Around MWC 297, we find an inner compact 
disk (r ou t~ 29 AU) and an outer ring at a distance of 200 AU. 
The 1 .3mm/2.7mm spectral index indicates that grain growth 
has proceeded in all disks in our sample. 

- A comparison between our data and previous results in TT 
and HAe stars show that the masses of the circumstellar disks 
around HBe stars are on average 5-10 times lower than those 
around Herbig Ae stars. We propose that disk photoevapo- 
ration is responsible for this behavior. In HBe stars, the UV 
radiation disperses the gas in the outer disk on a timescale of 
a few 10 5 yr. Once the outer part of the disk dissapears, the 
entire gaseous disk is photoevaporated on a short timescale 
(~10 5 yr), and only a dusty disk composed of large grains 
remains. 
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Table A.l. Table of fluxes for MWC 137. 






Wavelength (jim) 


rlux \jy) 


Flux error (ly) 


Beam size (") 


Reference / Flux type 


0.440 


0.019 


0.001 


0.8 


TYCHO-2 BAND.B 


0.550 


0.081 


0.001 


0.8 


TYCHO-2 BAND.V 


0.700 


0.017 


0.003 


0.9 


USNO-B BAND_RJOHNSON_MORGAN 


0.900 


0.012 


0.002 


0.9 


USNO-B BAND JJOHNSONJVIORGAN 


1.235 


0.499 


0.010 


2 


BANDJ.2MASS 


1.620 


0.428 


0.001 


1 


Morel 1978 BAND_HJOHNSON_MORGAN 


1.662 


0.749 


0.017 


2 


BAND.H_2MASS 


2.159 


1.500 


0.030 


2 


BAND.Ks.2MASS 


2.200 


1.719 


0.001 


1 


Morel 1978 BAND_KJOHNSON_MORGAN 


3.400 


3 





1 


Morel 1978 BAND_LJOHNSON_MORGAN 


8.280 


4.308 


0.01 


18.3 


BAND.AJVISX6C 


12.000 


10.9 


0.1 


30 


Hillenbrand 1992 (IRAS) 


12.130 


4.371 


0.010 


18.3 


BAND.CJVISX6C 


14.650 


2.713 


0.010 


18.3 


BAND.DJVISX6C 


21.340 


4.834 


0.010 


18.3 


BAND.EJVISX6C 


25.000 


25.7 


0.1 


30 


Hillenbrand 1992 (IRAS) 


60.000 


245.0 


1.0 


60 


Hillenbrand 1992 (IRAS) 


100 


519.0 


1.0 


120 


Hillenbrand 1992 (IRAS) 


350 


4 


1 


18.5 


Mannings 1994 


450 


1.66 


0.25 


18.3 


Mannings 1994 


750 


0.300 


0.060 


18.0 


Mannings 1994 


800 


0.370 


0.050 


17.0 


Mannings 1994 


1100 


0.095 


0.014 


18.8 


Mannings 1994 


1300 


0.090 


0.013 


19.8 


Mannings 1994 


1400 


0.0071 


6xl0~ 4 


3.05x1.78 


Our data 


2700 


0.0041 


2.0xl0- 4 


o.zxi.oi 


Our data 


7000 


2.35 10~ 3 


1.3xl0- 4 


2.21x1.66 


Our data 


13000 


0.0024 


5*10- 4 


5.18x3.19 


Our data 


36000 


5.7xl0- 4 


5xl0- 5 


0.3 


Skinner 1993 


60000 


<6.0*10- 4 




1.7 


Skinner 1993 


Table A.2. Extinction details for MWC 137. 


Wavelength 


Extinction 


Extinction from envelope model 


Observed flux Dereddened flux 


(jim) 


(mag) 


(mag) 




(mjy) (mjy) 


0.44 


7.80 


1.05 




19.2 25340.3 


0.55 


5.88 


0.79 




81.1 18294.1 


0.7 


4.42 


0.59 




17.3 1013.9 


0.9 


2.82 


0.38 




11.6 156.2 


1.235 


1.69 


0.23 




498.5 2372.9 


1.62 


1.09 


0.15 




428.1 1173.0 


1.662 


1.05 


0.14 




748.7 1969.6 


2.159 


0.69 


0.09 




1495.3 2821.1 


2.2 


0.67 


0.09 




1719.3 3182.7 
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Wavpl PTi dtfi ( nm\ 

Y Y Cl V l^l^llg Lll U-*'' l / 


Flux Clvt 

A ILIA \J J ) 


PIiiy prrnr f T\A 

A ILIA L11U1 \ J J J 


Rpfim si7P ( f, \ 


Reference / Flux type 




0.360 


0.043 


0.002 


1 


Reed 2005 BANDJJJOHNSONJVIORGAN 


0.440 


0.110 


0.001 


1 


AC2000.2 BANDJJJOHNSONJVIORGAN 


0.550 


0.166 


0.007 


1 


Hipparcos BAND.VJOHNSONJVIORGAN 


0.550 


0.219 


0.009 


1 


Corporon 1999 BAND JVJOHNSONJVIORGAN 


0.700 


0.208 


0.008 


1 


Zacharias 2005 BAND_RJOHNSON_MORGAN 


1.235 


0.580 


0.012 


2 


BANDJ.2MASS 




1.620 


0.328 


0.001 


1 


Morel 1978 BAND.HJOHNSON.MORGAN 


1.662 


0.762 


0.020 


2 


BANDJL2MASS 




2.159 


1.024 


0.020 


2 


BAND_Ks_2MASS 




2.200 


1.036 


0.001 


1 


Morel 1978 BAND_KJOHNSON_MORGAN 


3.400 


1.280 


0.001 


1 


Morel 1978 BAND.L JOHNSON JVIORGAN 


8.280 


3.261 


0.01 


18.3 


BAND.AJVISX6C 




12.000 


6.3 


0.1 


30 


Hillenbrand 1992 (IRAS) 




12.130 


4.2301 


0.01 


18.3 


BAND.CJVISX6C 




14.650 


1.84 


0.01 


18.3 


BANDJ3JVISX6C 




21.340 


2.15 


0.01 


18.3 


BANDJ2JVISX6C 




25.000 


8.2 


0.1 


30 


Hillenbrand 1992 (IRAS) 




60.000 


86.9 


0.1 


60 


Hillenbrand 1992 (IRAS) 




100 


315.0 


1.0 


120 


Hillenbrand 1992 (IRAS) 




450 


1.23 


0.22 


18.3 


Mannings 1994 




800 


0.142 


0.020 


17.0 


Mannings 1994 




1100 


0.058 


0.014 


18.8 


Mannings 1994 




1400 


<0.0010 




3.32x1.84 


Our data 




2700 


<3.0xl0- 4 




6.77x3.68 


Our data 




Table A.4. Extinction details for LKHa 215. 


Wavelength 


Extinction 


Extinction from envelope model 


Observed flux 


Dereddened flux 


(jim) 


(mag) 


(mag) 




(mJy) 


(mJy) 


0.36 


5.37 


0.78 




42.7 


6009.6 


0.44 


4.54 


0.66 




109.5 


7148.9 


0.55 


3.42 


0.50 




166.3 


3886.4 


0.55 


3.42 


0.50 




219.2 


5123.3 


0.70 


2.57 


0.37 




208.2 


2220.5 


1.235 


0.99 


0.14 




579.8 


1436.7 


1.62 


0.64 


0.09 




327.8 


589.0 


1.662 


0.61 


0.09 




761.9 


1337.3 


2.159 


0.40 


0.06 




1024.1 


1481.4 


2.20 


0.39 


0.06 




1036.0 


1482.1 
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Table A.5. Table of fluxes for R Mon. 



WflVPlPTlCTtfl ( nm\ 

Y Y Cl V l^l^llg Lll U-*'' l / 


Flux dvl 

A ILIA \ J V J 


PlllY PTTfir fTv^ 
A ILIA *J1 \J j } 


Beam size (") 


Rpfprenrp / FIiiy tvnp 


1.235 


0.213 


0.005 


2 




BANDJ.2MASS 


1.256 


0.027 


0.0010 


6 




Close 1997 


1.633 


0.25 


0.01 


6 




Close 1997 


1.662 


0.630 


0.015 


2 




BAND_H_2MASS 


2.159 


1.860 


0.040 


2 




BAND_Ks_2MASS 


2.210 


0.99 


0.01 


6 




Close 1997 


3.163 


11.17 


0.01 


6 




Close 1997 


3.600 


18 





12 




Audard 2007 


3.930 


19.69 


0.01 


6 




Close 1997 


4.290 


21.05 


0.01 


18.3 




BAND_B1_MSX6C 


4.350 


29.16 


0.01 


18.3 




BAND_B2_MSX6C 


4.674 


23.85 


0.01 


6 




Close 1997 


5.800 


30.41 


0.4 


12 




Audard 2007 


8.000 


31 





12 




Audard 2007 


8.280 


32.09 


0.01 


18.3 




BAND_A_MSX6C 


11.700 


46.59 


0.01 


6 




Close 1997 


12.000 


58.3 


0.1 


30 




Hillenbrand 1992 (IRAS) 


12.130 


47.18 


0.01 


18.3 




BAND_CJVISX6C 


14.650 


58.53 


0.01 


18.3 




BAND_D_MSX6C 


20.800 


103.1 


0.1 


6 




Close 1997 


21.340 


102.2 


0.01 


18.3 




BAND_E_MSX6C 


25.000 


139.0 


1.0 


30 




Hillenbrand 1992 (IRAS) 


52.560 


99 


4 


30 




Spitzer .tbl file 


69.660 


81 


3 


30 




Spitzer .tbl file 


93.600 


62.3 


2.0 


30 




Spitzer .tbl file 


350 


2 


0.300 


18.5 




Mannings 1994 


450 


1.56 


0.25 


18.3 




Mannings 1994 


800 


0.220 


0.040 


17.0 




Mannings 1994 


1 1U0 


O.U// 


0.018 


18.8 




Mannings 1994 


1200 


0.049 


10 


1 1 




Our data - 30m bolometer map 


1300 


0.079 


0.019 


19.8 




Mannings 1994 


1400 


0.0118 


0.0020 


0.72x0.33 




Our data 


2700 


0.0041 


5.0xl0- 4 


1.28x0.83 




Our data 


7000 


1. 26x10^ 


1.3xl0- 4 


2.21x1.66 




Our data 


13000 


7.8xl0- 4 


l.OxlO" 4 


5.18x3.19 




Our data 


20000 


5.9xl0- 4 


3.0xl0- 4 


1.4 




Skinner 1993 


36000 


4.0xl0- 4 


1.0xl0~ 4 


0.3 




Skinner 1993 


60000 


2.9xl0- 4 


l.OxlO- 4 


4.3 




Skinner 1993 


Table A.6. Extinction details for R Mon. 


Wavelength 


Extinction 


Extinction from envelope model 




Observed flux Dereddened flux 


(jim) 


(mag) 


(mag) 




(mJy) 


(mJy) 


1.235 


0.66 


1.36 




212.9 


390.8 


1.256 


0.64 


1.33 




27.0 


48.8 


1.633 


0.42 


0.87 




250.0 


368.3 


1.662 


0.41 


0.85 




629.7 


917.7 


2.159 


0.27 


0.55 




1863.5 


2386.2 


2.210 


0.26 


0.53 




990.0 


1256.2 


3.163 


0.15 


0.30 




11170.0 


12767.7 
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Table A.7. Table of fluxes for Z CMa. 
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Wavplpnotfi ( nm\ 

Y Y Cl V l^l^llg Lll U-*'' l / 


Flux dv^ 

L 1 UA \ J y J 


PIiiy error f T\A 

L 1 UA 1 KJL \ J ) 


Ream si7P (" \ 


Rpfprpnpp / Plnx tvnp 


0.440 


0.205 


0.001 


0.8 


TYCHO-2 BAND.B 


0.550 


0.474 


0.001 


0.8 


TYCHO-2 BAND.V 


0.700 


0.7 


0.3 


0.9 


USNO-B BAND_R_JOHNSON_MORGAN 


0.900 


0.9 


0.4 


0.9 


USNO-B BANDJJOHNSONJVIORGAN 


1.235 


4 





2 


BANDJ.2MASS 


1.620 


6 





1 


Morel 1978 BAND_HJOHNSON.MORGAN 


1.662 


8 





2 


BAND_H_2MASS 


2.159 


21 


3 


2 


BAND.Ks.2MASS 


2.200 


25 





1 


Morel 1978 BAND_KJOHNSON.MORGAN 


3.400 


52 





1 


Morel 1978 BAND_LJOHNSON_MORGAN 


3.800 


54.0 


4.0 


<0.1 


Koresko 1991 


4.290 


71.78 


0.01 


18.3 


BAND_B1_MSX6C 


4.350 


72.12 


0.01 


18.3 


BAND_B2_MSX6C 


4.800 


69.0 


6.0 


<0.1 


Koresko 1991 


8.280 


107.7 


0.1 


18.3 


BAND_A_MSX6C 


10.000 


109.0 


12.0 


<0.1 


Koresko 1991 


12.000 


126.6 


0.1 


30 


Oudmaijer 1992 (IRAS) 


12.130 


140.2 


0.1 


18.3 


BAND_C_MSX6C 


14.650 


159.1 


0.1 


18.3 


BAND_D_MSX6C 


20.000 


183.0 


5.0 


10 


ISO 


21.340 


202.8 


0.1 


18.3 


BAND_E_MSX6C 


25.000 


221.3 


0.1 


30 


Oudmaijer 1992 (IRAS) 


40.000 


265.0 


5.0 


20 


ISO 


50.000 


340.0 


10.0 


30 


Elia 2004 (ISO, chart estimate) 


60.000 


322.0 


0.1 


60 


Oudmaijer 1992 (IRAS) 


60.000 


365.0 


10.0 


40 


Elia 2004 (ISO, chart estimate) 


70.000 


390.0 


10.0 


40 


Elia 2004 (ISO, chart estimate) 


80.000 


410.0 


10.0 


60 


Elia 2004 (ISO, chart estimate) 


90.000 


420.0 


10.0 


80 


Elia 2004 (ISO, chart estimate) 


100 


354.0 


0.1 


120 


Oudmaijer 1992 (IRAS) 


100 


410.0 


10.0 


80 


Elia 2004 (ISO, chart estimate) 


120 


390.0 


10.0 


80 


Elia 2004 (ISO, chart estimate) 


140 


380.0 


10.0 


80 


Elia 2004 (ISO, chart estimate) 


160 


310.0 


10.0 


80 


Elia 2004 (ISO, chart estimate) 


180 


270.0 


10.0 


100 


Elia 2004 (ISO, chart estimate) 


350 


29 


1 


19.0 


Dent 1998 


450 


13.84 


0.2 


17.5 


Dent 1998 


450 


11.41 


0.1 


8.0 


Sandell 2001 


800 


1.964 


0.013 


16.0 


Dent 1998 


850 


1.5 


0.1 


15.0 


Sandell 2001 


1100 


0.710 


0.030 


18.7 


Dent 1998 


1300 


0.83 


0.01 


21.0 


Sandell 2001 


1400 


0.0260 


0.0010 


2.19x0.85 


Our data 


2700 


0.0085 


7.0xl0~ 4 


3.05x1.41 


Our data 


6917 


0.0021 


LOxlO" 4 


1.44x1.43 


Our data 


13350 


0.0021 


1.0xl0~ 4 


3.53x2.71 


Our data 


36000 


0.0022 


3.0X10- 4 


9.18x6.56 


Our data 



Table A.8. Extinction details for Z CMa. 



Wavelength 
ijim) 


Extinction 
(mag) 


Extinction from envelope model 
(mag) 


Observed flux 
(mJy) 


Dereddened flux 
(mJy) 


0.44 


5.09 


0.0 


204.6 


22253.1 


0.55 


3.84 


0.0 


474.0 


16279.5 


0.7 


2.88 


0.0 


689.6 


9818.2 


0.9 


1.84 


0.0 


890.4 


4850.9 


1.235 


1.11 


0.0 


3848.5 


10654.0 


1.62 


0.71 


0.0 


6075.0 


11728.6 


1.662 


0.69 


0.0 


8392.7 


15778.0 


2.159 


0.45 


0.0 


20774.5 


31437.6 


2.2 


0.44 


0.0 


24850.8 


37144.3 
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Table A.9. Table of fluxes for MWC 297. 



VV £l V GIGllgLll \JXfil) 


Flux (lv} 
i ilia x-f y j 


i iua ciivji \j y ) 




Reference / Flux type 




0.440 


0.012 


0.002 


0.9 


USNO-B BAND_BJOHNSON_MORGAN 


0.550 


0.046 


0.002 


<1 


Eisner 2004 BAND.V JOHNSON.MORGAN 


0.700 


0.300 


0.060 


0.9 


USNO-B BAND_RJOHNSONJVIORGAN 


0.791 


0.609 


0.005 


3 


BAND.LDENIS 




0.900 


1.040 


0.210 


0.9 


USNO-B BANDJJOHNSON.MORGAN 


1.228 


8 





3 


BANDJ.DENIS 




1.235 


6 





2 


BANDJ.2MASS 




1.620 


6 





1 


Morel 1978 BAND_H_JOHNSON_MORGAN 


1.662 


18 


2 


2 


BAND_H_2MASS 




2.159 


41 


4 


2 


BAND.Ks.2MASS 




2.200 


37 





1 


Morel 1978 BANDJCJOHNSONJVIORGAN 


3.400 


78 





1 


Morel 1978 BAND_L_JOHNS ON JVIORG AN 


4.290 


97.89 


0.01 


18.3 


oA.\ U_o l_lvloAOC 




4.350 


90.04 


0.01 


18.3 


DA.N l_-'_0 __1VT /\Uv 




8.280 


141.2 


0.01 


18.3 


RAMn A AA*\Yf\C 

DA.N L/_/\_IVlo AUC 




12.000 


159.0 


1.0 


30 


Hillenbrand 1992 (IRAS) 




12.130 


124.7 


0.01 


18.3 


BAND.CJVISX6C 




14.650 


104.8 


0.01 


18.3 


BAND.DJVISX6C 




21.340 


114.2 


0.01 


18.3 


BAND.EJVISX6C 




25.000 


224.0 


1.0 


30 


Hillenbrand 1992 (IRAS) 




60.000 


914.0 


1.0 


60 


Hillenbrand 1992 (IRAS) 




100 


1800.0 


10.0 


120 


Hillenbrand 1992 (IRAS) 




350 


4 





18.5 


Mannings 1994 




450 


2.46 


0.13 


18.3 


Mannings 1994 




600 


0.970 


0.110 


17.5 


Mannings 1994 




750 


0.840 


0.060 


18.0 


Mannings 1994 




800 


0.699 


0.013 


17.0 


Mannings 1994 




850 


0.681 


0.022 


18.0 


Mannings 1994 




1100 


0.452 


0.014 


18.8 


Mannings 1994 




1300 


0.432 


0.025 


19.8 


Mannings 1994 




1300 


0.3 


0.0 


3.14x3.02 


Manoj 2007 




1300 


0.175 


0.005 


1.1x0.43 


Our data 




2600 


0.149 


0.0050 


1.42x0.88 


Our data 




6917 


0.029 


0.0010 


1.95x1.64 


Our data 




13350 


0.0260 


2.0xl0~ 4 


3.92x3.06 


Our data 




36000 


0.01258 


2.0xl0~ 4 


7.2 


Skinner 1993 




60000 


0.00917 


2.0xl0- 4 


14.4 


Skinner 1993 




Table A.10. Extinction details for MWC 297. 


Wavelength 


Extinction 


Extinction from envelope model 


Observed flux 


Dereddened flux 


(jim) 


(mag) 


(mag) 




(mJy) 


(mJy) 


0.44 


9.24 






12.2 


60645.3 


0.55 


6.97 






45.8 


28036.0 


0.70 


5.23 






301.0 


37286.4 


0.791 


4.26 






608.9 


30758.9 


0.90 


3.34 






1041.4 


22567.3 


1.228 


2.02 






7726.1 


49860.6 


1.235 


2.01 






5645.3 


35818.6 


1.62 


1.30 






6420.2 


21181.2 


1.662 


1.24 






18009.4 


56620.4 


2.159 


0.82 






40507.0 


85900.8 


2.2 


0.79 






36926.7 


76571.6 
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Table A.ll. Table of fluxes for MWC 1080. 
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WavplpTidtfi (nm\ 


Flux f!v> 

± ILIA \J J ) 


PIiiy prrnr f T\A 

L 1 UA 1 I V / 1 \ J y J 


Ream <;i7p C r \ 


Reference / Flux type 




0.440 


0.022 


0.001 


0.8 


TYCHO-2 BAND.B 




0.550 


0.080 


0.003 


<1 


Eisner 2004 BAND_VJOHNSON_MORGAN 


0.550 


0.069 


0.001 


0.8 


TYCHO-2 BAND.V 




0.700 


0.120 


0.024 


0.9 


USNO-B BAND_RJOHNSONJVIORGAN 


0.900 


0.210 


0.040 


0.9 


USNO-B BANDJJOHNSON.MORGAN 


1.235 


1.650 


0.030 


2 


BANDJ.2MASS 




1.620 


1.921 


0.001 


1 


Morel 1978 BAND_H_JOHNS ON JVIORG AN 


1.662 


4 





2 


BAND_H_2MASS 




2.159 


8 





2 


BAND.Ks.2MASS 




2.200 


9 





1 


Morel 1978 BANDJCJOHNSONJVIORGAN 


3.400 


17 





1 


Morel 1978 BAND_L_JOHNSON_MORGAN 


4.350 


16.66 


0.01 


18.3 


BAND.B2.MSX6C 




8.280 


16.99 


0.01 


18.3 


BAND.AJVISX6C 




12.130 


17.91 


0.01 


18.3 


BAND.CJVISX6C 




14.650 


13.88 


0.01 


18.3 


BAND.DJVISX6C 




21.340 


16.12 


0.01 


18.3 


BAND.EJVISX6C 




60.000 


150.0 


1.0 


60 


Hillenbrand 1992 (IRAS) 




100 


246.2 


0.1 


120 


Hillenbrand 1992 (IRAS) 




350 


<24.55 


- 


18.5 


Mannings 1994 




450 


5 





18.3 


Mannings 1994 




600 


1.61 


0.22 


17.5 


Mannings 1994 




750 


0.774 


0.06 


18.0 


Mannings 1994 




800 


0.638 


0.014 


17.0 


Mannings 1994 




850 


0.660 


0.080 


18.0 


Mannings 1994 




1100 


0.250 


0.030 


18.8 


Mannings 1994 




1300 


0.237 


0.03 


19.8 


Mannings 1994 




1400 


0.0031 


2.0xl0~ 4 


1.28x0.79 


Our data 




2700 


<0.0017 




4.89x2.75 


Our data 




36000 


1.4xl0- 4 


2.0xl0- 5 


1.0 


Skinner 1993 




Table A.12. Extinction details for MWC 1080. 


Wavelength 


Extinction 


Extinction from envelope model 


Observed flux 


Dereddened flux 


(jim) 


(mag) 


(mag) 




(mJy) 


(mJy) 


0.44 


6.16 


4.24 




21.9 


6363.0 


0.55 


4.64 


3.19 




79.6 


5738.3 


0.55 


4.64 


3.19 




68.7 


4952.7 


0.70 


3.49 


2.40 




119.8 


2977.8 


0.90 


2.23 


1.53 




209.7 


1630.0 


1.235 


1.34 


0.92 




1652.3 


5662.8 


1.62 


0.86 


0.59 




1921.1 


4257.5 


1.662 


0.83 


0.57 




4152.4 


8911.5 


2.159 


0.54 


0.37 




7825.8 


12917.3 


2.20 


0.53 


0.36 




8696.4 


14141.4 
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Table A.13. Compilation of disks masses 1 . 



Name 


spectral 


Luminosity 


Teff 


Distance 


Star mass 


Age 


A 


Flux 


Disk mass 




type 


(L G ) 




(pc) 


/Ml \ 

(M Q ) 


(Myr) 


(mm) 


(mJy) 


/"A It \ 

(M ) 


AB Aur 


AO 


49 


9500 


144 


2.57 


3.87 


1.4 


85 


0.0278 


HD 100546 


B9 


45 


10500 


103 


2.53 


5.59 


1.3 


465 


0.0540 


HD 142527 


F6 


69 


6300 


200 


3.54 


1.11 


1.3 


1190 


1.4860 


HD 179218 /MWC 614 


B9 


500 


10000 


240 


5.02 


0.473 


1.3 


71 


0.0447 


HD 100453 


A9 


8 


7400 


114 


1.66 


14.1 


1.2 


270 


0.0709 


CD-36 10010B 


F8 


8 


6200 


84 


1.63 


9.49 


1.3 


142 


0.0345 


HD 139614 


A7 


11 


7800 


157 


1.83 


10.7 


1.3 


242 


0.1293 


HD 169142 


B9 


44 


10500 


145 


2.5 


5.69 


1.3 


169 


0.0389 


HD 142666 


A8 


8 


7500 


116 


1.65 


14.9 


1.2 


79 


0.0201 


HD 144432 


A9 


32 


7300 


200 


2.26 


3.87 


1.3 


37 


0.0369 


HD 150193 /MWC 863 


Al 


30 


9500 


150 


2.25 


6.10 


1.3 


45 


0.0143 


HD 163296 


Al 


36 


9500 


122 


2.49 


5.79 


1.3 


705 


0.1485 


HD 34282 


AO 


36 


9500 


400 


2.49 


5.79 


1.3 


110 


0.2259 


HD 35187 


A2 


34 


9100 


150 


2.43 


4.65 


1.2 


20 


0.0055 


UXOri 


A3 


42 


8600 


430 


2.43 


4.22 


1.2 


19.8 


0.0481 


CQTau 


F2 


5 


6800 


100 


1.49 


769 


1.3 


143 


0.0415 


MWC 758 


A3 


21 


8500 


200 


2.17 


6.26 


1.3 


72 


0.0471 


LkHo-259 


A9 


107 


7320 


850 


3.6 


1.17 


2.6 


<6 


<0.722 


WW Vul 


A3 


43 


8600 


550 


2.44 


4.19 


1.2 


9.1 


0.0362 


BF Ori 


A5 


34 


8300 


450 


2.34 


4.37 


1.3 


6 


0.0228 


SV Cep 


AO 


14 


8600 


715 


1.91 


11.7 


1.3 


7 


0.0459 


V376 Cas 


B5 


190 


13800 


600 


3.77 


2.66 


2.7 


<5 


<0.089 


LkHo-198 


B7 


210 


12000 


600 


3.51 


1.49 


2.7 


<5 


<0.1141 


VY Mon 


B8 


1400 


10700 


800 


6.9 


0.243 


1.3 


120 


0.6874 


HD 97048 


AO 


44 


10000 


150 


2.53 


4.61 


1.3 


452 


0.1305 


LkHo-234 


B5 


800 


14500 


1250 


5.01 


0.524 


1.3 


<20 


<0.181 


TOri 


A3 


83 


8500 


460 


2.92 


2.49 


1.3 


<6 


<0.021 


V380 Ori 


B9 


100 


9500 


460 


3.05 


2.00 


1.3 


<24 


<0.056 


LkHo-218 


B9 


107 


10700 


1150 


2.97 


2.77 


1.3 


<27 


<0.3906 


Elias 3-1 


A6 


21 


8000 


140 


2.06 


6.11 


2.7 


45 


0.1360 


BD+404124 


B2 


5900 


22000 


1000 


- 


- 


2.7 


<5 


<0.133 


VI 686 Cyg 


B5 


1820 


15500 


1000 


6.05 


0.312 


2.7 


<5 


<0.246 


BD+61154 


B8 


225 


11800 


650 


3.88 


1.43 


2.7 


11.2 


0.3509 


RRTau 


A4 


150 


8800 


800 


3.57 


1.34 


1.3 


<20 


<0.224 


HD 200775 


B2 


1900 


18500 


600 


6.9 


0.337 


1.3 


<6 


<0.007 


MWC 480 


A3 


25 


8500 


130 


2.19 


5.98 


1.4 


235 


0.0797 


I 1 U jt /UU 


VJW 


oc\ 
__w 


DUUU 


on 


Z.jO 


J.J7 


i . j>j 


1 1 7 

11./ 


U.UUj / 


MWC 297 


B1.5 


10700 


24000 


250 


9 


1 


1.3 


168.6 


0.0323 


ZCMa 


B8 


9200 


20000 


930 


12 


0.05 


1.4 


24.1 


0.2307 


R Mon 


B0 


5600 


25000 


800 


8 


1 


1.4 


8.3 


0.0093 


MWC 1080 


B0 


5600 


25000 


1000 


10 


1 


1.4 


2.1 


0.0036 


MWC 137 


B0 


15000 


26000 


1300 


14 


1 


1.4 


7.1 


<0.007 


VVSer 


AO 


50 


10000 


260 


2.68 


3.87 


1.3 


1.44 


0.0012 


LkHa 215 


B7.5 


5300 


22000 


800 


7 


0.4 


1.4 


<1.5 


<0.009 



1 The flux shown is the excess at 1.3 mm subtracting the free-free contribution in our sources. In MWC 297 we have used the flux of Mano i et al.l 
(2007). References for the mm fluxes are: Acke et al. 2004, Alonso-Albi et al. 2008, Di Francesco et al. 1997, Fuente et al. 2003, Fuente et al. 
2006, Henning et al. 1994, Henning et al. 1998, Mannings et al. 1997, Mannings et al. 2000, Natta et al. 1997, Natta et al. 2001, Pietu et al. 2003, 
Pietu et al. 2006, Sylvester et al. 2001. 



